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1. Introduction 

Relativistic fluid dynamics is an important model to understand various collective phenom- 
ena in astrophysics and heavy-ion collisions, although its theoretical foundation has not yet been 
established [[I]]. The relativistic Navier-Stokes theory is, for example, acausal and unstable and 
inadequate as the theory of relativistic fluids. The reason is that the irreversible currents (the shear 
stress tensor 7i^^ , the buUc viscous pressure n etc.) are linearly proportional to the thermody- 
namic forces (the shear tensor a^^, the expansion scalar Q etc.), with the proportionality constant 
named the shear viscosity coefficient T], the bulk viscosity coefficient etc. Thus, the forces have 
an instantaneous influence on the currents, which obviously violates causality and leads to insta- 
bilities. These problems aie solved by, for example, introducing retardation into the definitions 
of the irreversible currents, leading to equations of motion for these currents which thus become 
independent dynamical variables. The retardation effect is characterized by the relaxation time. 
Theories of this type are called causal dissipative relativistic fluid dynamics (CDR). In CDR, the 
irreversible currents and the thermodynamic forces aie no longer in a simple linear relation, and 
such fluids are called non-Newtonian. As a consequence, the transport coefficients for CDR can- 
not be computed with methods commonly used for Newtonian (Navier-Stokes) fluids, such as the 
Green-Kubo-Nakano (GKN) formula. 

Recently, a new microscopic formula to calculate the transport coefficients of CDR from time- 
correlation functions was proposed ||^ ^. This formula reproduces the ordinary results when it is 
applied to the classical Navier-Stokes theory and the diffusion equation. The consistency between 
this new formula and the results obtained from the Boltzmann equation was confirmed in Ref. 
[0, ^]. Since this formula is derived from quantum field theory, it will be applicable even to dense 
fluids, differently from the calculations based on the Boltzmann equation. 

The purpose of the present study is to calculate the transport coefficients of CDR with lat- 
tice QCD simulations by using the new formula. The calculations of the transport coefficients, 
in general, contain temporal-correlation functions which are very difficult to estimate in lattice 
simulations [^, ^. Thus, as a first attempt, we focus on a ratio between the shear viscosity and 
the corresponding relaxation time, T]/T;r, which is given only by static correlation functions. After 
defining the correlation functions between the energy-momentum tensor on the lattice, we calculate 
the ratio in quenched lattice simulations on 24^ x 4-16 lattices with jS^AT = 6.0, which corresponds 
to the temperature range 0.5 T/Jc ^ 1.8 where Tc is the critical temperature. 

This report is organized as follows: in section |2|, we introduce formulations of CDR and show 
that the ratio of transport coefficients can be expressed in terms of static correlation functions 
between the energy-momentum tensors. In section ^, the energy-momentum tensor is defined on 
the lattice by using a clover-shaped combination of gauge links. Results of lattice simulations are 
shown in section ^ and the summary is given in section |5[ 

2. Causal dissipative relativistic fluid dynamics 

We first choose gross variables which are necessary to extract the macroscopic motion of 
many-body systems. If the chosen variables are not enough, the derived fluid dynamics will show 
unphysical behaviors, such as instability and the divergent transport coefficients. 
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For ideal fluid, the energy-momentum tensor T^^ is a function only of the energy density £ and 
the fluid velocity , which is normalized asu^u^ = 1. Then, by applying a Lorentz transformation 
and using the definition of the energy density and pressure P, we obtain T^^ = {e +P)u^u^ —g^^P. 
Note that P is calculated by the equation of state. Since T^^ is conserved, we have 

d^T^'' = 0. (2.1) 

This is the relativistic Euler equation. 

For dissipative fluid, T^^ cannot be expressed only by e and u^. We represent this additional 
component by another second rank tensor n'^^. The most general T^^^ is, then, given by T^^ = 
{e +P)u^u^ — g^^ P + Yl^^ . Conventionally, n'^^ is expressed using the trace part IT and traceless 
part as Yl^^ = — {g^^ — u^u^)Yl. Finally T^^ is expressed as 

^ = (£ + P + U)u^'u'' - g^^P + n) + (2.2) 

and n and 7i^^ are the bulk viscous pressure and the shear stress tensor, respectively, satisfying the 
orthogonality condition u^Ti^^ = 0. In the traditional Landau-Lifshitz theory [||], the viscous terms 
are induced instantaneously by the corresponding thermodynamic force: 

n = -C0, 7r^^ = 2T]a^\ (2.3) 

where and rj are the bulk and shear viscosities, respectively. The thermodynamic forces 6 and 
a^^ are defined by 

d = d^u^, a'''' = ^(^d^u'' + d''u^-^{g^'''-u^u'')d^=A^'''^^d^us. (2.4) 

When we use these definitions of the viscous terms, we obtain the relativistic Navier-Stokes equa- 
tion. Because of the instantaneous production of the viscous terms, this equation contains sound 
propagations with infinite speed. 

In order to solve this problem, the retardation effect is taken into account by introducing relax- 
ation times T;r for the shear stress tensor and Tn for the bulk viscous pressure, respectively. Thus 
the viscous terms satisfying causality are given by 

Tnu^d^n + Tnnd+n=-^e, r^A^'''^^u"danxs + rKn^''e + n^''' = 2rio^'', (2.5) 

where Tn and are the relaxation times of IT and n^^, respectively. Here the projection operator 
^/ivA5 necessary to satisfy the orthogonality relation. These are the equations of CDR. One can 
easily check that the Navier-Stokes theory is reproduced in the vanishing relaxation time limit. 
The second terms on the l.h.s. come from the (de)compression of fluid cells which is important to 
implement stable numerical calculations with ultra-relativistic initial conditions [^. 

In fluid dynamics, transport coefficients are inputs which should be calculated from the un- 
derlying microscopic dynamics. As was discussed in the introduction, we cannot apply the GKN 
formula to CDR. The new formula is derived by using the projection operator method [^, 0]. The 
results are summarized as 

ri riGKN ^ 

j8(e-FP) ~/32/j3x(fov(x),fo-^'(0))' P 

C _ ^GKN 

j8(e + P) ~ j82/t/3x(fo^(x),fo-^(0))' J 



/32/j3x(f>'-^(x),f>'^(0))' ^ ■ 
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where " denotes operator, and we define IT = ^/^j T /3 — c^T and oA=A — Tr[pg^AJ with the 
equilibrium density matrix p^q. The inner product is defined by Kubo's canonical correlation, 

= ^Tr[p,,A(-a)B]. (2.8) 
Jo p 

Here '{]gkn and I^gkn are the shear and bulk viscosities of Newtonian fluids which are calculated 
using the GKN formula (or more precisely, using the Zubarev method). These quantities are given 
by the temporal (dynamical) correlation functions. 

One can see that the new transport coefficients are still calculated from the GKN formula with 
the normalization factors, which are, on the other hand, given by the static correlation functions. 
Thus, for example, the ratio of the shear viscosity and corresponding relaxation time is calculated 
only from the static correlation functions, 

r] _/J^^x(f--(x),f-«(0)) ^2.9) 



Tn{e+p) /J3x(r<'^(x),r<'^(0)) 



p 

e+P 



In the leading order of the weakly interacting bose gas, the above ratio becomes i^.p^ 
which becomes zero {\) at r = (r ^ oo) for massive bosons. In the following, we focus on this 
ratio and calculate it in quenched lattice simulations. 



3. Energy-momentum tensor on the lattice 



Let us consider the gluonic matter at finite T, and define the energy-momentum tensor for the 
SU(3) gauge theory in Euclidean space-time as. 



r^vW = 2tr [F^a{x)Fya{x)\ - ^5^v (l + ^) tr [Fpa[x)Fpa{x)\ , 



(3.1) 



where the trace is taken over color indices, and j3(g) is a beta function on the lattice [10]. In the 
standard approach, the field strength tensor squared on the lattice (without the summation over 
Lorentz indices) is defined from the Hermitian part of the plaquette as 



a^tr [F^vW^Mv W] + 0{a^) = 



1-^Retrf/^vW 



(3.2) 



where jS^^^^ = is a lattice gauge coupling. This is utilized to define e.g. the standard gauge 
action. However this does not tell us anything about the off-diagonal part of the energy-momentum 
tensor, T^v{p ^ v). Therefore, the following equality (valid only in the continuum theory with 
full 0(3) rotational symmetry) has been employed to calculate the correlations of the energy- 
momentum tensor: 



Tij{x)Tij{y)) = - [{T,{x)Tu{y)) - {Tn{x)Tjj{y))] , = 1,2,3). 



(3.3) 



It was however realized recently that this relation receives large errors due to lattice discretization 
[11]. Moreover, it does not give us a clue to calculate the correlation of 7]-4 (the denominator of the 
ratio in Eq. (2.9)) at finite T. 
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Alternative way to define the field strength would be to take the anti-Hermitian part of the 
plaquette, 



a^tr [F^,y{x)Fpa{x)\ +0{a^ 



24 



tr [Q^A^] - [Qpaix] - QlA^)] , (3.4) 



which can be used both for the first and second terms of the right hand side of Eq.( |3.1|) . Here we 
adopt a clover-shaped combination of the plaquette [12| 



1 



Q^iv{x) = - [[/^vW + t/v-Ai(x)+t/-^j-v(x) + t/-v/jW] 



(3.5) 



to respect the space-time symmetry. This definition naturally leads to (7)jv) = for /i / v. In our 
simulation, we use the energy-momentum tensor obtained from Eq.( |3.4[ ). 

In the Euclidean space-time, the Kubo's canonical correlation for the energy-momentum ten- 
sors appearing in Eq. ( ^^ ) becomes a susceptibility 



(3.6) 



where we have used the translation invariance both in spatial and temporal directions, and (• • •)r de- 
notes the thermal average at temperature T . With T = 1/ {aNt), V = {aNs)^ and / d^x f dx ^ <^^llx 
on the lattice, we can rewrite the static susceptibility G^y in the lattice unit with zero temperature 
subtraction as 



1 



1 



r=o. 



(3.7) 



where Nto means the temporal lattice size at T = 0. 



4. Results of lattice simulations 



We perform quenched lattice simulations employing a standard plaquette gauge action on a 
isotropic lattice of 24^ x Nt with Nt = 4— 16. The lattice coupling is taken to be jS^^^^ = 6.0, which 
corresponds to a = 0.093 fm with the Sommer scale ro = 0.5 fm [ 13]. The range of Nt corresponds 
to T/Tc ~ 0.5-1.8, where the critical temperature is located between Nt = 7 and N = 8. The zero- 
temperature subtraction is performed with N = 24. We generate pure gauge configurations by 
the pseudo-heat-bath algorithm and measure correlations using 1000-5000 configurations at every 
1000 trajectories after thermalization. Statistical errors are estimated by the jackknife analysis. 



In order to see the behavior of the energy-momentum tensor constructed from Eq. ( [3.4[ ), let us 
first show results of the trace anomaly, 



3P 



1 



T=0 



NM 



to 



X ^ 




(4.1) 



Figure [Jleft) shows temperature dependence of the trace anomaly together with the energy density 
and pressure calculated by the T-integral method [13]. Typical enhancement of {e — 3P) /T'^ around 
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Figure 1: Results of the trace anomaly, energy density and pressure (left) and the ratio of the susceptibilities 
Gxy{T)/GxA{T) (right) as a function of temperature. 



Tc, and the rapid (slow) increase of the energy density (pressure) can be seen. The off-diagonal parts 
of the energy-momentum tensor are found to be zero within the statistical error, {T^v)t — Oi (M / 
V). 



We define the averaged static susceptibilities Gxy and Gx4 from Eq.(|3.7|) as 



G^iT) = i (Gi2 + Gi3 + G23) , G,4(r) = i (Gi4 + G24 + G34) . (4.2) 

From the simulation, we found that both Gxy and Gx4 increase monotonically with temperature with 
similar values, so that the ratio Gxy/GxA shown in Fig. |I| (right) corresponding to r] / Xj^{e + P) is 
almost unity over the range of temperatures we have explored, 0.5 ^ T /Tc ^ 1.8. This behavior is 
in contrast to that expected from the weakly interacting bose gas mentioned at the end of sec.|2[ and 
is worth to be studied further. 



5. Summary 

We examined transport coefficients of causal dissipative relativistic fluid dynamics (CDR) in 
quenched lattice simulations. Based on the microscopic formulae proposed in Refs. [Q, ^], a ratio 
between the shear viscosity and the corresponding relaxation time, r\ / {Tji{e + P)) , was computed 
from the static correlation functions of the energy-momentum tensor. We calculated these static 
correlation functions in quenched lattice simulations on 24^ x 4-16 lattices with ^^^^ = 6.0, which 
correspond to the temperature range of 0.5 ^ T /T^- ^ 1.8. In this temperature region, the ratio stays 
constant and close to unity. 
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